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LONG  TERM  GOALS 

Our  long  term  goals  are  to  understand  the  dynamics  of  atmospheric  motions  on  scales  of  10 
kilometers  and  less  in  sufficient  detail  to  be  able  to  provide  a  consistent  subgrid  scale 
turbulent  closure  for  models  across  a  range  of  scales,  and  to  be  able  to  utilize  simulated 
variances  as  a  measure  of  forecast  predictability. 

OBJECTIVES 

The  objective  is  to  develop  a  consistent  closure  model  for  cloud  scale  turbulence  that  would 
allow  one  to  produce  a  simulation  at  any  desired  level  of  resolution,  with  the  results  of  lower 
resolution  simulations  being  approximately  similar  to  results  obtained  by  appropriate  spatial 
filtering  of  the  higher  resolution  simulation.  The  particular  objective  of  this  past  year  has 
been  to  make  progress  on  understanding  entrainment  at  the  top  of  cloudy  boundary  layers. 

APPROACH 

This  research  involves  the  utilization  of  the  high  resolution  turbulent  transport  codes 
developed  under  previous  ONR  support  to  explore  interesting  cloud  structure  features  which 
should  be  included  in  the  subgrid  parameterization  of  lower  resolution  models,  as  well  as 
progress  in  modeling  the  conservation  equations  governing  the  cloud  scale  variance  of 
temperature,  humidity,  and  velocity,  so  that  these  variables  may  be  included  in  larger  scale 
models. 

WORK  COMPLETED 

During  the  past  year  we  have  concentrated  on  understanding  the  results  of  our  large-eddy- 
simulation  (LES)  investigation  of  the  sensitivity  of  cloud-top  entrainment  to  a  number  of 
variables,  with  particular  emphasis  on  quasi-steady,  buoyantly-driven  boundary  layers. 
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Another  important  part  of  our  effort  over  the  current  year  has  been  continued  participation  in 
inter-model  comparisons  with  the  GEWEX  Cloud  Systems  Studies  boundary  layer  cloud 
modeling  working  group.  Our  LES  results  for  a  boundary  layer  with  shallow  cumulus  from 
BOMEX  were  included  in  the  1997  Seattle  workshop  (Siebesma,  1997). 

RESULTS 

Lewellen  and  Lewellen  (1997a,b)  use  large  eddy  simulations  to  strongly  support  the  rather 
unexpected  result  that  cloud  top  entrainment  is  often  controlled  by  large  eddy  transport 
rather  than  the  small  eddies  confined  to  the  local  mixing  zone.  In  simulations  of  quasi-steady, 
buoyantly-driven  boundary  layers,  the  turbulent  entrainment  rate  was  found  to  be  relatively 
insensitive  to  numerical  or  physical  alterations  which  affect  the  small  scale  mixing  dynamics 
near  the  inversion  but  leave  the  boundary  layer  scale  eddy  dynamics  largely  unchanged  (e.g., 
changes  in  grid  resolution,  subgrid  modeling,  or  modest  wind  shear).  The  entrainment  rate  is 
sensitive  to  changes  which  affect  the  large  scale  dissipation  dynamics  (e.g.,  low  Reynolds 
number)  as  one  might  expect  given  that  entrainment  and  dissipation  are  the  two  competing 
energy  sinks  in  these  systems. 

We  found  a  large  eddy  entrainment  efficiency,  T|e,  which  we  defined  as  the  fraction  of  the 
buoyant  energy  production  going  into  entrainment  rather  than  lost  to  turbulent  dissipation,  to 
be  fairly  constant  over  a  wide  range  of  quasi-steady,  buoyantly  driven  boundary  layers  with 
broadly  similar  large  eddy  dynamics.  As  we  have  implemented  it,  T|e  is  similar  (but  not 
identical)  to  the  flux  partitioning  measure  of  Stage  and  Businger  (1981).  Since  this  measure 
of  entrainment  cost  is  a  global  one,  its  implementation  requires  knowledge  of  boundary  layer 
properties  such  as  cloud  fraction  and  depth;  on  the  other  hand  it  allows  us  to  avoid  defining  a 
local  buoyancy  jump  at  Zj,  which  requires  making  assumptions  about  mixing  at  cloud  top. 
This  efficiency  increases  from  approximately  0.3  for  a  surface  driven  dry  boundary  layer  with 
a  stable  top,  to  -0.35  for  a  top  driven  dry  boundary  layer  with  a  stable  top,  and  up  to  0.5  for  a 
quasi-steady  boundary  layer  with  an  unstable  cloud  top  with  the  cloud  fraction  reduced 
below  one. 

Entrainment  results  as  parameterized  by  our  entrainment  efficiency,  T|e,  and  the  more 
traditional  A  parameter  (The  ratio  of  the  entrainment  velocity  over  the  characteristic 
convective  velocity,  w*,  times  the  Richardson  number  based  on  the  virtual  potential 

temperature  jump  (A0V)  across  the  inversion,  the  mixed  layer  thickness  (zi),  and  w *2  in  the 
layer  (i.e.,  Bretherton  et  al.  1997))  are  summarized  for  many  simulations  in  Fig.  1.  Where 
possible,  results  from  two  different  one  hour  averaging  periods  are  included  to  give  some  idea 
of  the  time  variability  within  a  single  simulation.  The  entrainment  efficiency  is  seen  to  be 
fairly  constant  while  the  A  parameter  varies  by  more  than  an  order  of  magnitude.  There  is  in 
general  a  negative  contribution  to  the  buoyancy  flux  throughout  the  boundary  layer  arising 
from  entrainment,  as  the  warm  air  from  above  the  inversion  is  mixed  throughout  the 
boundary  layer  down  to  the  ground.  Despite  differences  in  inversion  jumps  and  flux  profiles, 
the  boundary  layer  dynamics  for  the  cases  considered  are  broadly  similar,  being  dominated 


by  large  eddies  spanning  the  boundary  layer.  The  relative  constancy  of  follows  as  long 
as 


A 


Fig.  1.  Summary  of  simulation  results  showing  entrainment  efficiency  r|£  versus  the  more 
traditional  A  parameter  as  defined  m  the  text.  The  simulation  cases  are  labelled  as  follows: 
smoke  cloud  with  surface  heating  (•);  smoke  cloud  with  cloud-top  radiational  cooling  with 
different  optical  depths  (♦);  radiatively  cooled  smoke  cloud  with  cooling  height  varied  (o); 
radiatively  cooled  smoke  cloud  with  varying  wind  shear  (□);  radiatively  cooled  water  cloud 
with  varying  humidity  jump  across  the  inversion  (A);  radiatively  cooled  water  cloud  with 
L/cp  varied  (A);  radiatively  cooled  water  cloud  with  varying  surface  heat  flux  (0).  Those 
points  overscored  with  "x"  indicate  unstable  cloud  tops. 

the  large  scale  mixing  involved  in  the  entrainment  process  is  the  dominant  limiting  factor 
setting  the  entrainment  rate. 


IMPACT 


A  consistent  quantitative  model  of  cloud  top  entrainment  is  important  to  any  model  which 
involves  cloud  dynamics.  In  addition  to  the  navy's  operational  forecasting  interest  in  clouds, 
an  understanding  of  cloud  dynamics  on  this  scale  is  also  a  central  issue  in  modeling  global 
climate  change. 

TRANSITIONS 

We  sent  information  to  Steve  Burk,  Marine  Meteorology,  Naval  Research  Laboratory, 
responding  to  his  request  for  possible  parameterization  of  the  inversion  thickness  which  might 
be  used  in  radar  refractivity  estimates. 

RELATED  PROJECTS 

The  LES  code  developed  under  ONR  support  has  been  modified  and  used  to  model  aircraft 
contrails  for  NASA  (Le  we  lien  and  Lewellen,  1996),  and  to  model  the  turbulent  interaction  of 
a  tornado  with  the  surface  for  NSF  (Lewellen  et  al.  1997). 
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